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We calculate the ground state properties of fermionic dipolar atoms/molecules in a ID double- 
tube potential by using the Luttinger liquid theory and the Density Matrix Renormalization Group 
calculation. When the external field is applied near a magic angle with respect to the double-tube 
plane, the long-ranged dipolar interaction can generate a spontaneous correlation between fermions 
in different tubes, even when the bare inter-tube tunneling rate is negligibly small. Such interaction- 
induced correlation strongly enhances the contrast of the interference fringes, and therefore can be 
easily observed in the standard time-of-flight experiment. 



I. INTRODUCTION 



Recent progress of ultracold atoms has made it possi- 
ble to study strongly correlated physics in a much wider 
parameter range. One of the most important subjects 
is the one-dimensional (ID) physics, and many interest- 
ing phenomena, including Tonks-Girardeau gas [l], Lut- 
tinger liquid (LL) behavior ^], and polaronic effects in 
Bose- Fermi mixture [3]], etc. have been theoretically pro- 
posed or experimentally observed. However, due to the 
short-range nature of atomic interaction, it is not easy to 
study how the interaction between particles in different 
ID tubes can bring different many-body effects. In tra- 
ditional solid state systems, on the other hand, the long- 
ranged Coulomb interaction between electrons does in- 
duce several important many-body features in low dimen- 
sional multi-component systems: for example. Coulomb 
drag between ID double wires [4|, inter- wire/ well coher- 
ence [5i,6|], and spontaneous ferromagnetism (or exciton 
condensation) in double layer quantum Hall systems Q- 
Following these extensive studies in the the solid state 
systems, it is therefore very interesting to investigate 
how the long-ranged dipolar interaction between ultra- 
cold atoms/molecules [8|9| can bring emergent many- 
body physics in a spatially separated multi-component 
system, where the first-order and second-order correla- 
tion function can be observed easily in a time-of-flight 
(TOF) measurement. We note that several pioneering 
works have been reported to explore the exotic quan- 
tum phases of bosonic dipolar atoms/molecules [10] in 
the multi-layer/tube systems, but results for fermionic 
dipoles seem not well-explored yet to the best of our 
knowledge. 

In this paper, we investigate a ID double-tube system 
(see Fig . [Ija), and the similar setup for bosonic atoms in 
Ref. [Ill), where fermionic uoI'ay molecules 9} (or mag- 
netic dipolar atoms, ^^Cr Q) are loaded in the ground 
state of the confinement potential. When no external 
field is applied, fermions arc nonintcracting and have no 
correlation between fermions in the two tubes. Therefore 
the time-of-fiight (TOF) image (imaging light is along the 
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FIG. 1: (Color online) (a) The double-tube system considered 
in the paper. Tilted arrows indicate the dipole moment of 
dipolar atoms/molecules. An external field (E for the electric 
field or B for the magnetic field) is applied with an angle 9 
in the x — y plane. The big leftward arrows represent the 
imaging light for taking a TOF image in the fast expanding 
(i.e. y — z) plane, (b) and (c) are the calculated momentum 
distribution {N±(ky,kz), proportional to the TOF image), 
based on the DMRG result of the lattice model (see Eq. (I14|l 
and (jlSp below), (b) is the noninteracting result, and (c) is 
the interacting result with 6 = 9c. All the parameters used 
here are described in the text. The dashed and solid curves in 
(d) are the amplitude profile along the horizontal line {ky — 0) 
of (b) and (c) respectively. 



tube direction) should be structureless (see Fig. [TJb)). 
However, when an external electric/magnetic field is ap- 
plied and tilted near a magic angle, 9c — cos~^ a/1/3, the 
intra-tube dipolar interaction is strongly reduced, leav- 
ing only a repulsive and long-ranged inter-tube interac- 
tion. We show that such exotic interaction can induce 
the first- order correlation between fermions in the two 
tubes, which can be easily measured from the interfer- 
ence fringes (see Fig. [He)). We note that unlike the 
interference between bosons |lllll3 |. the interference be- 
tween such spatially separated fermionic particles has no 
counterpart in classical waves, and is induced purely by 



the strongly correlated effect. 

To investigate this problem systematically, in the fol- 
lowing we will first apply the celebrated LL theory with 
the renormalization group (RG) study in the weakly in- 
teracting limit, and then use Density Matrix Renormal- 
ization Group (DMRG) method in the strong interac- 
tion regime. The weak or strong interaction regime is 
defined by the ratio between dipolar interaction to the 
kinetic energy. Results from both calculations assure the 
existence of such interaction-induced inter-tube correla- 
tion, and provide a rich quantum phase diagram. For 
the convenience of later discussion, we will introduce the 
pseudo-spin up/down index for the upper/lower tube, so 
that the conventional definition of magnetic orders can 
be used to identify the quantum phases in the present 
system. For example, the ground state with spontaneous 
inter-tube correlation can be understood as an in-plane 
ferromagnetic order, {tpl{x)4'\{x)) j^ 0, where ipy^{x) is 
the fermionic operator in the upper/lower tube. As a 
result, the inter-tube correlation can be understood as a 
ID planar ferromagnetic (FM) state with the pseudo-spin 
polarized in the x — y plane of spin space. This is exactly 
the ID version of the well-known spontaneous ferromag- 
netism in double layer quantum Hall systems l7|| . Similar 
problem about ID ferromagnetism of itinerant fermions 
is also an important subject discussed in the literature 
|13l |. Note that the pseudo-spin defined here has nothing 
to do with the original spin of fermionic dipoles. 



II. HAMILTONIAN AND LUTTINGER LIQUID 

THEORY IN THE WEAKLY INTERACTING 

LIMIT 

Throughout this paper, we will assume that parti- 
cles are loaded in the lowest subband of each tube 
with a transverse confinement wavefunction (j)s{y,z) = 
^i^e"'^ +(z-sd) )/2R g^j^i^ ^ Gaussian radius R (here 

s = ±i=t/J,is the pseudo-spin index). The result- 
ing effective ID system Hamiltonian then can be written 
tohe H = Hq + Hi, where Hq is the kinetic energy: 



Hq 



E 



dxi}l{x)-^dl'!j)^{x) 



(1) 



with m being the mass of dipolar particles, and Hi is the 
interaction energy: 
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(2) 



Here V\s-s'\ix) is the dipolar interaction between 
molecules in the same tubes (Vb) and/or in dif- 
ferent tubes (Vi), obtained by integrating out 
the transverse degree of freedom: V|5_s/| (x) = 
J dri^± f dr2,± \(f>s{ri^±)\'^ |0s'(r2a)|^ Vd (ri -r2). 



where Vrf(r) = _D^(1 — 3(r ■ £')^)/|rp is the bare dipolar 
interaction with D being the electric dipole moment in 
c.g.s unit. E is the unit vector along the external electric 
field. Since in general the electric dipole interaction is 
much stronger (and tunable) than the magnetic dipole 
interaction, in the the rest of this paper we will use 
polar molecules as the underlying particles for further 
discussion. Extension to the magnetic dipolar atoms is 
straightforward. 

Defining Vo/i(fc) = J dxVo/i{x) e~^''^ to be the Fourier 
transform of interaction, we can calculate their zero mo- 
mentum (fc = 0) value analytically by integrating over 
the transverse confinement wavefunction, (j)s{rj_): 
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(4) 



It is very easy to see that when 6 ^ 6c = cos^^ ^^1/3 the 
intra-tube interaction is almost zero while the inter-tube 
interaction is still finite and positive. Such interesting 
kind of interaction matrix element cannot be realized in 
the traditional solid state system and therefore may bring 
some physics not predicted before. 

In the weakly interacting limit, we can apply the 
standard LL theory [ij] by linearizing the band energy 
around the two Fermi points, i/cf, and dividing fermions 
into the left /right movers (i.e. tjjs{x) = J2r=±''i'r,s{x), 
where r = ± is the chiral index). The resulting low 
energy effective Hamiltonian can be rewritten to be, 
Hcs = i?LL + Hi, where 



i?LL = VF 



/ dx'Y^ V)J ,, (x) f -^rh- kp j i'r,s (x) 

^ r,s ^ ^ 

1 ■ / f^ 

2 E ^^sfs' / dxpr^s{x)pr',s'{x) 



(5) 



is the LL Hamiltonian with the density operator, 
Pr,s{x) = ipl ^{x)t[^r,s{x) and the forward scattering am- 
plitude: gl'l, = Vls-s'\{0) - 9\\Sr-r'Ss,s'; 



Hi = g± 






dx4i*j^,,{x)4'^ 



(x)V'L,_,(x)V'+,,(x) (6) 



is the backward scattering between particles of different 
tubes and different chiralties. Here vp is the Fermi ve- 
locity, and g\\n = VQ/i{2kp) is the backward scattering 
amplitude. According to the standard Luttinger theory 
[lj|, the linearized band structure about the two Fermi 
points is justified only when the interaction strength is 
much smaller than the Fermi energy, i.e. when 



i?3 



< 



2?Ti 



(7) 



It is well-known that the LL Hamiltonian, i?LLi can 
be diagonalized exactly via a Bogoliubov transformation. 



while the backward scattering term, Hi , cannot be diag- 
onaUzed in the same basis. As a result, we have to use 
the standard one-loop renormalization group (RG) calcu- 
lation [14| to investigate when the later term is relevant 
in the low energy limit, and how it renormalizes the for- 
mer one (-//ll) in different parameter regime. The bare 
Luttinger exponents are given by the initial system pa- 
rameters. 



Kp/a = 



TTVp 



law 



TTVF + Vp/o 



bw 



(8) 




where Vp/^ = Vb(0) ± ^1(0). Within the one-loop RG, 
only Ka and the backward scattering, g±, are renormal- 
ized [ijl , following the RG equations below: 



dl 
where 
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9± \ dg± 

7ru„ / at 



^{VF + K/vr)2 - (K/7r)2 



K.). (9) 



(10) 



is the collective mode velocity of the (pseudo-)spin-mode 
sector. / = In A is the scaling parameter with A being the 
shortest length scale (or the largest momentum scale) in 
the present system. For a given bare system parameter, 
{K„,g±), Eq. © then shows how they can flow to a fixed 
point, {K^,g*^), which determines the true low energy 
physical properties of the double-tube system. 



III. PHASE DIAGRAM IN THE LUTTINGER 
LIQUID THEORY 

To obtain a quantum phase diagram in the LL the- 
ory, we have to use the RG result to calculate the scaling 
exponent, a, of various correlation function in the long- 
distance limit, i.e. (O^(x)O(O)} ~ a;~^+" as |a;| -+ 00, 
where 0{x) is the order parameter operator of inter- 
est. For a given system parameter, the dominant or- 
der parameter is then determined by the most slowly 
decaying correlation function (or the largest and posi- 
tive a), known as the quasi-long-ranged order in ID sys- 
tems. Similar to the standard Luttinger liquid theory for 
other ID systems [IJj, we investigate the following can- 
didates of order parameters lj|, because scaling expo- 
nents for other more complicated kinds of order param- 
eters are always smaller than the ones below and there- 
fore becomes negligible. The definitions of order param- 
eters we considered in this paper includes (the associated 
scaling exponent is shown after the definition): Wigncr 

crystal (Owe = Es V'+,>+,-s'^-,-s^-,s' "wc = 2- 
^Kp), charge density wave (Ocdw = Z^s '^+,sV'-,s: 
acDW = 2 — Kp — Ka), axial (pseudo-)spin density 



wave (OsDW, = Y.s,s' ^l,s'i'+,ai'-,s'^ asDW, = 2- 
Kp — Kcr), planar (pseudo-)spin density wave (Osdw^,^ = 



FIG. 2: (Color online) Quantum phase diagrams in terms 
of the tube separation relative to the tube radius, d/R, and 
the tilted angle, O/tt, in the LL theory. {a,)kpR = 0.1 
(filled circle), (b)fcF-R = 0.2 (filled square) and kpR = 0.4 
(open circle). The dashed line indicates the magic angle, 
6c = cos~^ vl/3 « 0.3047r. (c)-(e) are cartoons to represent 
the axial SDW, polarized triplet superfluid, and the planar 
SDW phases respectively as defined in the text (also con- 
sistent with the literature [1J|). The filled and open circles 
here indicate the particle and hole fluctuation respectively. 
The the horizontal ellipses in (d) indicate pairing between 
fermionic particles, while the vertical ellipses in (e) indicates 
an inter-tube correlation (i.e. there is an uncertainty for the 
particle position between the two tubes) . The average particle 
density distribution is uniform for SDWa;,y phase (e), while 
it is staggered for SDWz phase (c), which has no inter-tube 
correlation. 



= 2 - i^„ - A'; 



and 



here o"^,'^'^ is Pauli matrix element), singlet superfluid 
{Oss = Es'SV'+,>- -s, ass = 2 - Kp^ - K„), unpo- 

larized triplet superfluid {OtSo = J2s '4'+,s4' s: Q^tSo — 

2 — K^^ — Ka), and polarized triplet superfluid {Ots2s = 

^J+Jj_^s for s - ±1/2, ars,^ ^2-K-^-K-^). In the 
low energy limit, the exponent a for each phase is then 
obtained from the renormalized Luttinger exponents, AT* 
(obtained via solving Eq. ©) and AT* — Kp. Note that 
K* is not renormalized within the one-loop RG scheme 
and therefore is the same as the initial(bare) value ob- 
tained from the system parameters. 

After calculating the scaling exponents of various or- 
der parameters mentioned above, we find that only three 
of them can be the dominant phases in the system we 
consider here. They are axial (pseudo-)spin density wave 
(SDW2), planar (pseudo-)spin density wave (SDWj;_j^), 
and polarized triplet superfluid (TS±). The charge den- 
sity wave (CDW) has the same scaling exponent as 
SDW2, but it can be shown that it still decay fast than 
the later case by a logarithmic correction after integrat- 
ing along the RG flow [1^]. In this paper, we will still 
consider both of them to be the dominant phase, but only 
discuss the physics related to the SDW2 order for simplic- 
ity. From the definition of scaling exponent above, it is 
easy to see that the phase boundary between the planar 
(pseudo-)spin density wave (SDWj;.y), and the polarized 



triplet superfluid (TS±) is determined by K* = Kp = 1, 
while the phase boundary between the axial (pseudo- 
)spin density wave (SDWz) and the planar (pseudo-)spin 
density wave (SDW^^^j^) is by K* = 1. From the expres- 
sion of the bare LL exponent in Eq. ([8]), we can find the 
first phase boundary is then determined by the following 
simple equation: 



Vp^Vo + Vi 



Voi2kf 



(11) 



However, when considering the phase boundary be- 
tween SDWz and ST)Wx,y, we have to calculate the renor- 
malized Luttinger exponent, K*, from the RG equation 
in Eq. ([9]) . Integrating Eq. ^ directly, it is easy to show 
that the solution connecting to K* = 1 and g^ = (re- 
quired as a fixed point of Eq. ©, dK^/dl = dgj_/dl = 0) 
is given by 



K„ 8 \TTVa 



(12) 



To obtain the functional dependence of the interaction 
strength on this RG line, we can insert the bare form of 
K„ in Eq. ([8]) and take the weak interaction limit. The 
leading order terms then give the following condition for 
the phase boundary between SDWz and SDW^j^j,: 

= 5± + .911 - Wkp) + Vi{2kF). (13) 

Since each term in Eqs. pT|) and p^ is proportional 
to D^ , the phase boundary lines determined above are 
independent of the interaction strength (more precisely, 
true only in the weak interaction limit). As a result, 
the phase diagram of LL theory can be completely deter- 
mined by the following three dimensionless parameters: 
d/i?, 0, and kpR. 

In Figs. [2Ia) and (b), we show the calculated phase 
diagram as a function of d/R and 0, for three differ- 
ent values of particle density, kpR. For a given inter- 
tube separation, d, the system is dominated by the ax- 
ial (pseudo-)spin density wave (SDWz) in the regime of 
large tilted angle 6. Such state can be also understood as 
a "phase-locked" state with staggered modulation of the 
particle density (see Fig. [2Ic)), similar to the ID elec- 
tron gas in the semi-conductor double wire system p^ . 
This phase results from the fact that when is large, 
the intra-tube interaction is repulsive and stronger than 
the inter-tube interaction, making a crystalized density 
distribution in each tube. At the same time, the back- 
ward scattering, (7_l, scales to be infinite (i.e. g^ ^ oo 
because K* < 1, see Eq. ([Q])), and therefore opens a 
single particle excitation gap. Such divergence of 5j_ also 
leads to a vanishing iiT* according to Eq. ^ . This is why 
the scaling exponent of SDWz becomes much larger than 
SDWa;_j, in the low energy limit (see above for the defini- 
tion of scaling exponents for each phase) . In the regime 
of small 6, on the other hand, the intra-tube interaction 
becomes attractive, leading to a triplet superfluid {TS±) 



(Fig. Hid)). It is easy to see from the order parameters of 
the above two phases, there is no inter-tube correlation, 
and therefore no visible interference fringes in the fast 
expanding plane, similar to the nonintcracting result in 
Fig. mh). 

However, when 6 is in the intermediate range and close 
to Oc ^ 0.3047r, the ground state becomes the planar 
(pseudo-)spin density wave (SDW^^^y). In this param- 
eter regime, the intra-tube interaction becomes smaller 
than the inter-tube interaction (i.e. |Vo(0)| < Vi(0), see 
Eqs. ©-(HI)), and therefore the backward scattering be- 
comes irrelevant in the low energy limit (5^ — > because 
K* > 1, see Eq. ©), leading to a uniform density dis- 
tribution along the tube. Besides, the order parameter, 
OsDW:r : has implied that at each position on the tube, 
fermionic particles has equal and non-zero probability to 
be found in the upper and the lower tube (or say, there 
is an uncertainty for a particle to be in one tube or the 
other, although the average probability in these two tubes 
are the same), showing an interaction induced inter-tube 
correlation (with a periodic modulation of the relative 
gauge phase, see Fig. (He)). Although within the LL 
theory, such inter-tube correlation is a quasi-long-ranged 
order, we still expect to observe it in the interference 
pattern (similar to Fig. [Ijc)) in a finite size system after 
integrating along the tube direction, if only the scaling 
exponent, a, is not too small. For typical parameters of 
polar molecules, say OH molecules (to ^ 17 a.m.u. and 
the largest dipole moment is D ~ 1.68 Debye), we can 
consider R ~ 0.1 /xm, d ~ 0.5 //m, and kp ^ ir x 10 
fim^^, and hance the obtained asDW^ can be of the 
order of 0.1 or more, showing that the probability to ob- 
serve such inter-tube correlation is not small. We note 
that this planar (pseudo-)SDWa;^y phase with a sponta- 
neous inter-tube correlation is a completely new phase, 
and has no counterpart even in the traditional solid state 
(semi-conductor or metal) double-wire system. 



IV. FERROMAGNETISM IN STRONGLY 
INTERACTING REGIME 

In the strongly interacting regime (i.e. the interaction 
energy is of the same order of or even larger than Fermi 
energy), LL theory fails to give correct low energy physics 
and hence we apply DMRG method to numerically study 
the ground state properties. For the convenience of nu- 
merical calculation, from now on we will consider 6 = 9c 
only to eliminate the intra-tube interaction completely. 
Besides, we assume an optical lattice is applied along the 
tube direction (x) so that the system Hamiltonian can 
be written in a single band lattice model: 



i?iat = -tYl (4«S+i,« + '^•^•) " '^51^^ 
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],B 
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+U^nj^^nj^l + - Y^ hj^. 



S"']2, — S 



(14) 
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where clj^s and fij^s 



^t 



«i,sS,« 



here are the fermionic 
field operator and density operator of the site j and 
the pseudo-spin s =t / [. Here i, J, C/, and V are 
the inter-site tunneling in the same tube, the inter- 
tube tunneling, the inter-tube onsite repulsion, and 
the inter-tube nearest-neighboring-site repulsion respec- 
tively. Note that we introduce inter-tube tunneling J 
here just for the completeness of numerical calculation, 
but we will just concentrate on the physics in the limit of 
zero J after the calculation, i.e. in the limit of zero inter- 
tube tunneling. The numerical values of these parame- 
ters can be easily calculated by integrating over the on- 
site Wannier function, and can be tuned in a wide range 
by changing the external field strength, lattice strength, 
and/or lattice spacing. Since it is not our purpose in this 
paper to provide a numerical comparison with any cer- 
tain experimental setup (there is unfortunately no such 
system available yet), we will not calculate their absolute 
values in this paper, but will present all of our results in 
dimensionless parameters instead. For simplicity, we also 
have neglected the interaction of longer (next-nearest- 
neighboring sites) range. 

In Fig. [3] we first show the ground state energy, Eq, 
as a function of V , by taking J = 0, U/t = 3 and 
the total particles, N — 10, in a ID lattice of length 
L = 20. Note that by taking J = 0, particle num- 
bers in each tube is a conserved quantity and there- 
fore we can label the system eigenstates by (7V| , ^i ) ■ 
There are three different regimes of interest: in regime I 
{V < 1.317), Eg increases linearly as a function of V, and 
the numbers of particles in the two tubes are the same 
(iV| — Ni — 5) in the ground state. Detailed analysis 
shows that the ground state waevfunction is the same 
as the SDW^^^y phase in the weak interaction limit dis- 
cussed above. In regime III {V > 1.321), Eg becomes 
independent of V, because all particles are moved to ei- 
ther one of the tubes to minimize the strong inter-tube 
interaction energy. This phase can be described as an 
axial ferromagnetic state by breaking the Z2 (spatial in- 
version) symmetry of the Hamiltonian, Eq. (|14[) . Inside 
the narrow window of regime II (1.317 < V < 1.321), 
the ground state energy changes step by step as the par- 
ticle number difference changes, indicating a series of 
first order phase transitions. In the thermodynamic limit 
{N,L ^ 1), it is reasonable to expect that the ground 
state energy should become a smooth curve in regime II, 
i.e. an infinitesmall inter-tube tunneling (J — > 0) can 
easily mix many degenerate states and lead to a coher- 
ent ground state, |\l/c), which has finite single particle 
inter-tube correlation between fermions in the two spa- 
tially separated tubes. Using the language of magnetism, 
we can define the pseudo-spin magnetization density to 

he M = jf J2j J2s.s' ^s,s'a}j,sO'j.s'^ ^^"^ ™ch spontaneous 
inter-tube correlation can be understood as an in-plane 
ferro-magnetization, i.e. (^cl-'^xl^c) 7^ 0. 

To confirm this conjecture, we further investigate re- 
sults with a finite but small inter-tube tunneling, J, 
which is equivalent to an effective " magnetic field" along 
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FIG. 3: (Color online) Ground state energy of Eq. (fTl)) (for 
J = 0) as a function of V/t for different particle numbers in 
the two tubes, (A''|, Ny). Regime I, II, and III are the SDwa;,y, 
canted ferromagnetic (or inter-tube correlated) phase, and ax- 
ial ferromagnetic phase respectively. The detailed definition is 
in the text. Inset shows the magnetic susceptibility x (defined 
in the text) as a function of inter-tube hopping, J. 



the X direction in the pseudo-spin space. In other words, 
we can study how the system magnetization responds to 
such small perturbation by calculating the "magnetic" 
susceptibility, x — d {"^ g\J^ x\^ c) / d J . As shown in the 
inset of Fig. [3l x tias a tendency to diverge as J ^ 
when V is tuned to be inside the regime II, strongly in- 
dicating an in-plane FM order. Combining the fact that 
the average numbers of particles in the two tubes are 
also different in this regime (i.e. {^c\-M.z\^g) 7^ 0, as 
indicated from Fig. [3]), the ground state in regime II is 
best understood as a pseudo-spin canted ferromagnetism, 
i.e. pseudo-spin is uniformly polarized with a tilted an- 
gle about the pseudo-spin z axis. This is another striking 
new phase predicted in this paper, breaking a continuous 
C/(l) symmetry in ID system. Different from the pla- 
nar SDW discussed in the weakly interacting limit, the 
pseudo-spin FM state obtained here in the strongly in- 
teracting regime has a uniform (rather than periodic os- 
cillation) magnetization. From the experimental point of 
view, this order parameter indicates a strong interference 
fringes as shown in Fig. [Ijc). Unfortunately, so far we 
have not been able to calculate systems of larger size with 
an efficient code, and therefore the demonstration of the 
existence of such symmetry breaking state in the ther- 
modynamic limit may need more justification. However, 
to the best of our knowledge, our system is probably 
the first physical candidate (rather than a mathemati- 
cal model) to realize the ferromagnetism of ID itinerant 
fermions: it is not be restricted by the celebrated Lieb- 
Mattis theorem [16| or Mermin-Wigner theorem [17| due 
to the nontrivial long-ranged and anisotropic dipolar in- 
teraction. Numerical results for larger system size and 
more complete phase diagram of this double-tube sys- 
tem should be an interesting direction to explore in the 
future. 



V. INTERFERENCE PATTERN IN THE 
TIME-OF-FLIGHT EXPERIMENT 

Since the theme of this paper is to demonstrate the 
existence of an interaction-induced inter-tube correla- 
tion between ID dipolar fermions (true for both weak 
and strong interaction regime), the best experimental 
evidence is from the interference pattern between these 
fermions. Theoretically, we can calculate the momen- 
tum distribution function of the ground state and then 
integrate along the tube direction to obtain the TOF im- 
age in the fast expanding (jj — z) plane, -/Vj^(kj^) (here 
k_L = {ky,kz), see Fig. [T]). Since particles can fly fast 
from the tube center after the strong confinement poten- 
tial is released, such TOF image should eventually be- 
come equivalent to momentum distribution after a long 
time of flight. After integrating out the contribution 
along the tube direction (x) , we can obtain the following 
expression for A'j^ (kj^ ) : 



A^j_(k_L) = AnR^Ne 



-|ki|"«" 



(l + 2(X,)cos(fc,d)H5) 



where we have used a Gaussian type-wavefunction (with 
radius R) to approximate the initial confinement wave- 
function, (l}sir±), and (Mx) = (^c\Mx\^c) is the inter- 
tube correlation. When the inter-tube tunneling, J, is 
small and the interaction is weak, the calculated inter- 
tube correlation {{A4x)) is so small that the obtained 
momentum distribution (i.e. TOF image) is close to a 
broad Gaussian-like function without any structure (see 
Fig. [Ub) and the dashed line in Fig. [^d)). On the other 
hand, when the direction and the strength of dipolar in- 
teraction is tuned to the regime II of Fig. [31 the long- 
ranged inter-tube interaction strongly enhances the inter- 
tube correlation {{Mx)) and hence the contrast of inter- 
ference fringes. In Fig. [TJc) (and the solid line in Fig. 
[TJd)), we show a typical calculated interference pattern 
(i.e. A^j^(k)) with parameters in this regime: U/t — 3, 
V/t = 1.318, and R/d = 0.167 for quarter-filling (L = 20 
and N = 10). We find that even the same single particle 
inter-tube tunneling rate {J/t = 0.1) is used in the cal- 
culation of both Fig. [Tfb) and (c), the obtained interfer- 
ence patterns are totally different: the later has strong 
contrast in the interference fringes due to the interac- 
tion effect. Unlike the usual interference fringes between 
bosonic particles, the fermionic dipoles here do not have 
a condensate in each tube, and therefore has no counter- 
part in any classical waves. 



As for the quasi-long-ranged order predicted by the 
Luttinger liquid theory in the weakly interacting regime, 
one should measure the TOF image in the x — z plane by 
using the imaging light perpendicular to the double tube 
plane (y). As a result both SDW^ and SDW^^^y phases 
can show some kinds of Bragg peaks along the x-axis due 
to the periodic distribution of density and/or phase fiuc- 
tuations, while the TS±i phase does not have such fea- 
ture. Furthermore, SDWz phase has a single particle gap, 
while SDWa;_y and TS±i do not have such gap. This gap 
can be easily measured by light scattering spectroscopy 
and hence becomes a way to distinguish these two density 
wave phases. We note that although the inter-tube corre- 
lation of SDWx,y phase is periodically oscillating and the 
correlation function decays as a power-law in large dis- 
tance, the finite tube length can still make it possible to 
measure some residual inter-tube correlation, which does 
not exist in the other two phases. According to the re- 
sults of bosonic atoms 11,18], the noise correlation of the 
TOF image for different lengths of tube can also provide 
the experimental value of the Luttinger exponent. 



VI. SUMMARY 

In this paper, we find an interaction-induced inter-tube 
correlation between spatially-separated fermionic polar 
molecules (or dipolar atoms) in the ID double-tube po- 
tential. Such correlation does not exist if the interaction 
between fermions is too weak, showing a clearly many- 
body effect in a ID system. We use both Luttinger liq- 
uid theory (with proper renormalization group method) 
to study the phase diagram in weak interaction limit, 
and use Density Matrix Renormalization Group Method 
to calculate the order parameters in strong interaction 
regime. This new phenomenon predicted in this paper 
can be observed in the first order interference pattern, 
and has no counterpart either in the classical waves or in 
the ID solid state systems. 
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